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Abstract. We study graphical mean curvature flow of complete solutions de- 
fined on subsets of Euclidean space. Wc obtain smooth long time existence. 
The projections of the evolving graphs also solve mean curvature flow. Hence 
this approach allows to smoothly flow through singularities by studying graph- 
ical mean curvature flow with one additional dimension. 
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O 1. Introduction 

O 

\Q Results. We start by stating a simplified version of our main result, which holds 

for bounded domains. Let us consider mean curvature flow for graphs defined on a 
relatively open set 



o 

(N 



(1.1) ft = (J fit x {t} c x [0, oo). 

i>0 

Then we have 



Theorem 1.1 (Existence on bounded domains). Let A C R n+1 be a bounded open 
set and u a : A — > R a locally Lipschitz continuous function with u a (x) — > oo for 
x — > x € dA. 

Then there exists (fi,u), where fi C R n+1 x [0, oo) is relatively open, such that 
u solves graphical mean curvature flow 

Dz 



u = 7W • div \-j===\ 4 nfi\(fi x{0}). 
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The function u is smooth for t > and continuous up to t = 0, Oo = A, u(-, 0) = uq 
in A and u(x,t) — > oo as (a;, t) — > 90, w/iere 90 is the relative boundary of m 
x [0,oo). 

Such smooth solutions yield weak solutions to mean curvature flow. To describe 
the relation, we use the measure theoretic boundary <9 M O t as introduced in Sec- 
tion A. We have the following informal version of our main theorem concerning the 
level set flow: 

Theorem 1.2 (Weak flow). Let (A,uo) and (Cl,u) be as in Theorem 1.1. Assume 
that the level set evolution o/<9O does not fatten. Then it coincides with (<9 M O t ) t > . 

For the general version of our existence theorem see Theorem 8.2. Theorem 9.1 
is our main result concerning the connection between the smooth graphical flow 
and the weak flow (in the level set sense) of the projections. In general, we do not 
know whether the solutions (O, u) are level set solutions. We notice, however, that 
such a statement would imply uniqueness of (ft, u) in Theorem 8.2. 

The previous theorems also provide a way to obtain a weak evolution of a set 
E C R n+1 with E = dA for some open set A: Consider a function uq : A — > M 
as described in Theorem 8.2, for example uq(x) := dist ^ + \x\ 2 , and apply our 
existence theorem. Then we define as the weak evolution of E the family (dttt)t>o 
with the notation from above. 

Illustrations. We illustrate our main theorems by some figures. In the description, 
we assume for the sake of simplicity that Q t — E t . 



In Figure 1 we study the evolution of a graph over B x (0) (drawn with thick lines), 
that is asymptotic to the cylinder S™ x R (drawn with grey lines) . The thinner lines 
indicate how the graph looks at some later time. We remark that it continues to be 
asymptotic to the evolving cylinder, which collapses in finite time. As we prove in 
Theorem 8.2, the evolving graph does not become singular and it has to disappear 
to infinity at or before the time the cylinder collapses. Theorem 9.1 implies that 
the evolving graph and the evolving cylinder disappear at the same time. Notice 
that near the singular time, the lowest point moves arbitrarily large distances in 
arbitrarily small time intervals. 

Figure 2 illustrates a graph over a set that develops a "neck-pinch" at t = T. 
This is projected to lower dimensions. For t /*T, the graph splits above the "neck- 
pinch" into two disconnected components without becoming singular. The thinner 
lines illustrate the graph for t > T. The rest of the evolution is similar to the 
situation above. 




FIGURE 1 . Graph over a ball 
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Figure 2. Graph over a set that develops a "neck-pinch" 




Figure 3. Graph defined initially over an annulus 

Next, we consider a rotationally symmetric graph over an annulus, centered at 
the origin, see Figure 3. The inner boundary of the annulus converges to a point 
as t /• T. At t = T a "cap at infinity" is being added to the evolving graph. This 
cap moves down very quickly. By comparison with compact solutions we see that 
u(0, t) is finite for any t > T. This is illustrated with thin lines. Finally, once again 
the evolution becomes similar to the evolution in Figure 1. 




Figure 4. Domain with nontrivial topology 
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Similarly, when a graph over a domain as in Figure 4 evolves, "caps at infinity" 
are being added at the times when the small "holes" shrink to points. 

Strategy of proof. In order to prove existence of smooth solutions, we start by 
deriving a priori estimates. The proof of these a priori estimates is based on the 
observation that powers of the height function can be used to localize derivative 
estimates in space. Then the result follows by applying these estimates to approx- 
imate solutions and employing an Arzela-Ascoli-type theorem to pass to a limit. 

The connection between singularity resolving and weak solutions is obtained 
as follows: We observe that the cylinder (dflt x K) t acts as an outer barrier for 
graphu(-, t). Furthermore, since graphu(-, t) — R converges to the cylinder as R — > 
oo, we conclude that graph u(-,t) does not detach from the evolving cylinder near 
infinity. 

Literature. The existence of entire graphs evolving by mean curvature flow was 
proved by K. Ecker and G. Huiskcn [11] for Lipschitz continuous initial data and 
by J. Clutterbuck [6], T. Colding and W. Minicozzi [8] for continuous initial data. 
K. Ecker, G. Huisken [10] and N. Stavrou [29] have studied convergence to homo- 
thetically expanding solutions, J. Clutterbuck, O. Schniircr, F. Schulze [5] and A. 
Hammerschmidt [20] have investigated stability of entire solutions. 

Many authors have worked on weak formulations for mean curvature flow, e. g. 
K. Brakke [3], K. Ecker [9], L. C. Evans, J. Spruck [12, 13, 14, 15], Y. Chen, Y. 
Giga, S. Goto [4] and T. Ilmanen [25]. In what follows we will refer as weak flow 
to level set solutions to mean curvature flow in the sense of Appendix A, see also 
[4, 12, 21]. 

Smooth solutions and one additional dimension have been used by S. Altschuler, 
M. Grayson [1] for curves to extend the evolution past singularities and by T. 
Ilmanen [24] for the e-regularization of mean curvature flow. 

Several people have studied mean curvature flow after the first singularity. We 
mention a few papers addressing this issue: J. Head [21] and J. Lauer [26] have 
shown that an appropriate limit of mean curvature flows with surgery (see G. 
Huisken and C. Sinestrari [22] for the definition of mean curvature flow with surgery) 
converges to a weak solution. T. Colding and W. Minicozzi [7] consider generic 
initial data that develop only singularities that look spherical or cylindrical. In 
the rotationally symmetric case, Y. Giga, Y. Seki and N. Umeda consider mean 
curvature flow that changes topology at infinity [17, 18]. 

The height function has been used before in [19] to localize a priori estimates for 
Mongc- Ampere equations. 

Organization of the paper. The classical formulation X = —Hv of mean curva- 
ture flow does not allow for changes in the topology of the evolving hypersurfaces. 
Hence in Section 2 we introduce a notion of graphical mean curvature flow that 
allows for changing domains of definition for the graph function and hence also 
changes in the topology of the evolving submanifold. 

We fix our geometric notation in Section 3 and state evolution equations of 
geometric quantities in Section 4. 

The key ingredients for proving smooth existence are the a priori estimates in 
Section 5 that use the height function in order to localize the estimates. 

In Section 8 we prove existence of smooth solutions. That result follows from 
combining the Holder estimates of Section 6 and the compactness result that we 
prove in Section 7 (a version of the Theorem of Arzcla-Ascoli). In Section 9 we 
discuss the relationship of our solution and the level set flow solution; we prove 
Theorem 9.1. Finally, we include an appendix that summarizes some of the results 
used in Section 9. 
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Open problems. Wc wish to mention a few open problems: 

(1) What is a good description of solutions disappearing at infinity? 

(2) If the projected solutions or a connected component of the complement 
become symmetric, e. g. spherical, does the graph pick up that symmetry? 

(3) What are optimal a priori estimates? 

(4) Is the solution (fl, u) unique? 

(5) Does the level set solution of graph u fatten? Is this fattening related to 
that of the level set solution of dAl 

Acknowledgment. We want to thank many colleagues for their interest in our 
work and inspiring discussions: G. Bellettini, K. Ecker, G. Huisken, T. Ilmanen, H. 
Koch, J. Metzger, F. Schulze, J. Spruck and B. White. Some of these discussions 
were possible due to invitations to Barcelona, Berlin, Oberwolfach and Potsdam. 

2. Definition of a solution 

Definition 2.1. 

(i) Domain of definition: Let fl C R n+1 x [0, oo) be a (relatively) open set. Set 
Q t := 7r R „+i (fl n (R n+1 x {t})), where 7r R „+i : E"+ 2 -> R n+1 is the projection 
to the first n + 1 components. Notice here that the first n + 1 components of 
the domain fl are spatial, while the last component can be understood as the 
time component. 

Observe that for each fixed t the section fit C R n+1 is relatively open, 
(ii) The solution: A function u: fl — > R is called a classical solution to graphical 
mean curvature flow in fl with continuous initial value u - f!o^ R> */ 

ue ^:(n\(flox{o}))ncL(S]) 

where we recall the definition of the spaces below and 

{ u = JTT\DuP ■ div ( . Du ) fflft\m x{0}), 
(MCF) < Vll VVl + l^l 2 / 

\u(-, 0) = u in flo- 

(Hi) Maximality condition: A function u: fl — > R fulfills the maximality con- 
dition if u > — c for some c e R and if w|n n (K"+ 1 x[o.T]) « s proper for every 
T>0. 

An initial value u ■ O — > R, flo C is said to fulfill the maximal- 

ity condition if w: flo x [0, oo) — > M defined by w(x,t) := Uo(x) fulfills the 
maximality condition, 
(iv) Singularity resolving solution: A function u: fl — > R is called a singular- 
ity resolving solution to mean curvature flow in dimension n with initial value 
u : fl -> R if 

a) fl and flo are as in (i), 

b) u is a classical solution to graphical mean curvature flow with initial value 
u as in (ii) and 

c) u fulfills the maximality condition. 

(v) We do not only call u a singularity resolving solution but also the pair (fl,u) 
and the family (M t ) t > with M t = graph u(-,t) C R n+2 . 

Remark 2.2. 

(i) Note that the domain of definition will depend on the solution. 

The dimensions seem to be artificially increased by one. This is due to 
the fact that we wish to study the evolution of (dflt)t>o, which in the smooth 
case, see Remark 9.9 (v) ; is a family of n- dimensional hypersurfaces in R n+1 
solving mean curvature flow. 
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(ii) If£l = K n+1 then condition (ii) in Definition 2.1 coincides with the definition 
in [11]. 

We avoid writing a solution as a family of embeddings X : M — > K™ +2 as 
in general, the topology of M is not fixed when fl t becomes singular. 

We expect similar results for other normal velocities, for example, if u is a 
singularity resolving solution for the normal velocity Sk in dimension n then 

u = a/1 + \Du\ 2 ■ S k [u] «Sl\(f! x{0}), 

where Sk[u] denotes the k-th elementary symmetric function of the n + l 
principal curvatures of graph u(-,t) C R n+2 and ft is as in Definition 2.1 (i). 
(Hi) a) The maximality condition implies that u tends to infinity if we approach 
a point in the relative boundary dfl. It also ensures that u(x,t) tends to 
infinity as \x\ tends to infinity. 

Hence the maximality allows us to use the height function u for localizing 
our a priori estimates. 

b) Our maximality condition implies that each graph 

M t = graph u{-,t) C R n+2 

is a complete submanifold. 

c) If u fulfills the maximality condition then u (x) := u(x,0) also fulfills the 
maximality condition. 

d) The maximality condition prevents solutions from stopping or starting sud- 
denly. Furthermore, in general restrictions of the domain £1 of a singularity 
resolving solution (Ct, u) do not provide other singularity resolving solu- 
tions, i. e. for general open sets B C W l+1 x [0, oo), the pair (ft n B,u\b) 
does not fulfill the maximality condition. 

(iv) It suffices to study classical solutions to mean curvature flow to obtain singu- 
larity resolving solutions. Nevertheless, this allows to obtain weak solutions 
starting with dQo by considering the projections of the evolving graphs. 

3. Differential geometry of submanifolds 

We use X = X{x, t) = (^ a )i< Q < n+2 to denote the time-dependent embedding 
vector of a manifold M" +1 into R™ +2 and j- t X = X for its total time derivative. 
Set M t := X(M, t) C W l+2 . We will often identify an embedded manifold with 
its image. We will assume that X is smooth. Assume furthermore that M n+1 is 
smooth, orientable, complete and dM n+1 = 0. We also use that notation if we have 
that situation only locally, e. g. when the topology changes at spatial infinity. 

We choose v = v{x) = (^ Q ) 1<Q ,<„ +2 to be the downward pointing unit normal 
vector to M t at x. The embedding X(-, t) induces at each point of M t a metric 
(9ij)i<i. j<n+i and a second fundamental form (/ijjOixi.jXn+i- Let (<7 lJ ) denote 
the inverse of {g%j). These tensors are symmetric and the principal curvatures 
(Aj)i<i< n +i are the eigenvalues of the second fundamental form with respect to 
that metric. As usual, eigenvalues are listed according to their multiplicity. 

Latin indices range from 1 to n + 1 and refer to geometric quantities on the 
surface, Greek indices range from 1 to n + 2 and refer to components in the ambient 
space M n+2 . In ]R n+2 , we will always choose Euclidean coordinates with fixed e n+2 - 
axis. We use the Einstein summation convention for repeated upper and lower 
indices. Latin indices are raised and lowered with respect to the induced metric or 
its inverse (g^), while for Greek indices we use the flat metric {g a p)i< a ,[3<n+2 = 

(Ml<a,/3<n+2 ofM"+ 2 . 

Denoting by (•, •) the Euclidean scalar product in we have 

9ij = X,j) = X^SafiX^j, 
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where we use indices preceded by commas to denote partial derivatives. We write 
indices preceded by semi-colons, e.g. hij-k or v-k, to indicate covariant differentia- 
tion with respect to the induced metric. Later, we will also drop the semi-colons 
and commas, if the meaning is clear from the context. We set X." = X" and 

;ij ,ij ij ,k' 

where 

Itf - y^U.L; I .'/,/.. 9ij,l) 
are the Christoffel symbols of the metric (gij). So becomes a tensor. 

The Gaufi formula relates covariant derivatives of the position vector to the 
second fundamental form and the normal vector 

(3.2) = h tjl y. 

The Weingarten equation allows to compute derivatives of the normal vector 

(3.3) i£ = h\X« k . 

We can use the Gaufi formula (3.2) or the Weingarten equation (3.3) to compute 
the second fundamental form. 

Symmetric functions of the principal curvatures are well-defined, we will use the 
mean curvature H = Ai + . . . + A n+1 and the square of the norm of the second 
fundamental form \A\ 2 = \\ + . . . + A^ +1 . 

Our sign conventions imply that H > for the graph of a strictly convex function. 

The space 

C k,a:k/2, a /2 denotes the 

space of functions for which up to fc-th 
derivatives are continuous, where time derivatives count twice, these derivatives are 
Holder continuous with exponent a in space and a/2 in time and the corresponding 
Holder norm is finite. The space Ci oc (fl) consists of the functions u: fl — > K which 
are in C k (K) for every K g O. We use similar definitions for other (Holder) spaces. 
Finally, we use c to denote universal, estimated constants. 

4. Evolution equations for mean curvature flow 

Definition 4.1. If M is given as an embedding and a graph, we use 77 = (0, ...,0,1) 

to denote the vector e n +2- The definitions of v, H and \A\ 2 are as introduced in 
the previous section. We denote the induced connection by V and the associated 
Laplace- Beltrami operator by A. 

We define v = (—ri a v a ) 1 and u = rj a X a . The function u can be regarded as 
a function defined on a subset of W l+l x [0, 00) or as a function defined on the 
evolving manifold M . It should be clear from the context which definition of u is 
being used. 

Theorem 4.2. Let X be a solution to mean curvature flow. Then we have the 
following evolution equations. 

(|-A)«=0, 

(i-A)v=-v\A\ 2 -l\Vvf, 
(£-A)\A\ 2 = -2|VA| 2 + 2|A| 4 , 
(f t - A) \\7 m A\ 2 < - 2 |V m+1 A| 2 

+ c(m, n) ■ \ V ™ A \ ■ I VM I ' I V3A \ ■ I yfej4 l » 

i-\-j-\-k—m 

(£ t -A)g<- 2kg 2 - 2^v- 3 (\7v, vg), 

where g = (f\A\ 2 = x \ v i |A| 2 and k > is chosen so that kv 2 < \ in the domain 
considered. 
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We remark that whenever we use evolution equations from this theorem, we 
consider u as a function denned on the evolving manifold. 

Proof. See [9, 11]. □ 

5. A PRIORI ESTIMATES 

The following assumption shall guarantee that we can prove local a priori es- 
timates for the part of graph u where u < 0. Notice that, via considering the 
evolution given by u — a (where a is a constant abbreviating the Spanish word 
"altura"), this is equivalent to obtain bounds in the set where u < a. 

In this section we will consider the set Cl = {u < 0}. More precisely, we will 
work under the following assumption: 

Assumption 5.1. Let Cl C R n+1 x [0, oo) be an open set. Let u: Cl — > R be a 

smooth graphical solution to 

u = y/1 + \Du\ 2 ■ div [ — 1 = Du | in VIC) (R n+1 x (0,oo)) . 

W 1 + \ Du \ 2 J 

Suppose that u(x,t) — > as (x,t) — > (x ,t ) € dCl. Assume that all derivatives of 
u are uniformly bounded and can be extended continuously across the boundary for 
all domains Cl (~l x [0, T]) and that these sets are bounded for any T > 0. 

Remark 5.2. 

(i) Assumption 5.1 is fulfilled for smooth entire solutions u to graphical mean 
curvature flow that fulfill u > L > 1 outside a compact set when we restrict u 
toCl = { (x, t) e R n+1 x [0, oo) : u(x, t) < 0} . 

(ii) The approximate solutions R in Lemma 8.1 fulfill Assumption 5.1 for L > 
0. 

(Hi) The following a priori estimates extend to the situation when 

ft = {(x, t) : u(x, t) < a} 

for any n€l instead o/0. We only have to replace u by (u — a) below, e. g. 
in Theorem 5.3. 

(iv) The boundedness assumption of the sets follows from the properness of the 
function u. 

Theorem 5.3 (C 1 -estimates). Let u be as in Assumption 5.1. Then 

2 ^ 2 

vu < max uu 

t=0 
{i»<0} 

at points where u < 0. 

Here and in what follows, it is often possible to increase the exponent of u. 
Proof. According to Theorem 4.2, w := vu 2 fulfills 

tit — vu 2 + 2vuu, 
Wi =ViU 2 + 2vuui, 

Wij =Viju 2 + 2vuuij + 2vUiUj + 2u(viUj + VjUi), 
( f t - A) w =u 2 ( f t - A) v - 2v\Vu\ 2 - 4u(Vw, Vu) 

= u 2 {-v\A\ 2 - 2 -\Vv\ 2 ) - 2v\Vu\ 2 - 4 (^V«, VvVu 

< -u 2 v\A\ 2 < 0. 

The estimate follows from the maximum principle applied to w in the domain where 
u < 0. □ 
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Remark 5.4. If the reader prefers to consider a positive cut-off function (—it), we 
recommend to rewrite Theorem 5.3 as an estimate for v ■ (—u) 2 . 

Corollary 5.5. Let u be as in Assumption 5.1. Then 

v < max vu 2 

t=a 
{«<o} 

at points where u < — 1. 

Remark 5.6. Similar corollaries also hold for higher derivatives. We do not write 
them down explicitly. 

Remark 5.7. For later use, we estimate derivatives of u and v, 
\Vu\ 2 ^aX^Xfnp = v* (s aP - v a v !i ) Vf3 = l-v- 2 <l 
and, according to (3.3), 

\\7v\ 2 = ((-ffei/T 1 )/'' ((V)" 1 ). = v^aX^g^Xfnp < v*\A\ 2 

<v 2 <p\a\ 2 = v 2 q. 

So we get 

\{Vu,Vv)\ < \Vu\ ■ \ Vv\ < v 2 \A\ < vVQ. 

Theorem 5.8 (C 2 -estimates). Let u be as in Assumption 5.1. 

(i) Then there exist A > 0, c > and k > (the constant in ip and implicitly in 
Q), depending on the C 1 -estimates, such that 

tu A g + Xu 2 v 2 < sup Xu 2 v 2 + ct 

t = 
{«<0} 

at points where u < and < t < 1 . 
(ii) Moreover, if u is in C 2 initially, we get C 2 -estimates up to t = 0: Then there 
exists c > 0, depending only on the C 1 -estimates, such that 

u 4 G < sup u A g + ct 

t=0 
{«<0} 

at points where u < 0. 

Proof. In order to prove both parts simultaneously, we underline terms and factors 
that can be dropped everywhere. We get the first part if we consider the underlined 
terms and the second part if we drop those and set A = 0. 
We set 

w := tu 4 g + Xu 2 v 2 

and obtain 

to = y^g_ + Atu 3 gu + tu 4 g + 1\v 2 uu + 2Xu 2 vi), 
Wi = 4:tu 3 gui + tu^gi + 2Xv 2 uui + 2Xu 2 Wi, 
Wij = 4:tu 3 guij + tu A gij + 2Xv 2 uuij + 2Xu 2 Wij + 12tu 2 guiUj 

+ 4tu 3 (giUj + gjUi) + 2Xv 2 UiUj + 2Xu 2 ViVj + 4Xvu(uiVj + UjVi), 

tu 3 \/g = -Vw - 4:tu 2 g\7u - 2Xv 2 \7u - 2Xuv\7v, 
u 
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(| - A) w < y^g_ + tu* (-2kG 2 - 2^v- 3 (\7v, Vg)) + 2Xu 2 v (-v\A\ 2 - ||Vw| 2 ) 

- 12tu 2 g\Vu\ 2 - 8tu 3 (Vg,Vu) - 2Av 2 |Vu| 2 - 2Xu 2 \Vv\ 2 

- 8Xuv(Vu,Vv). 

In the following, we will use the notation (Vro, b) for generic gradient terms for the 
test function w. The constants c are allowed to depend on sup{|ti| : u < 0} (which 
does not exceed its initial value) and the C 1 -estimates which are uniform as we 
may consider v ■ (u — l) 2 in Theorem 5.3. In case (i), it may also depend on an 
upper bound for t, but we assume that < t < 1. That is, we suppress dependence 
on already estimated quantities. 

We estimate the terms involving Vg separately. Let e > be a constant. We 
fix its value blow. Using Remark 5.7 for estimating terms, we get 

-2^ii 4 v~ 3 (Vv, VG) = -2^" / Vu, -Vw - 4iu 2 gVu-2Xv 2 Vu-2XuvVv) 

v 6 \ u I 

< (Vw,b)+8t^-g\A\+4X^\u\v\A\+4^-\Vv\ 2 

< CVw, b) + etu 4 g 2 + eXu 2 v 2 \A\ 2 + Xu 2 \Vv\ 2 ■ 4^ + c(e, A), 

v 

-8iw 3 (V£, Vu) = - 8 ( Vu, 1 Vic - Atu 2 gVu - 2Xv 2 Vu - 2XuvVv 
\ u 

< (Vw, b) + 32t_u 2 g + WXv 2 + 16X\u\v 3 \A\ 

< (Vw, b) + et_u A g 2 + eXu 2 v 2 \A\ 2 + c(e, A). 

We obtain 

(| - A) w < u A g + t_u A g 2 (-2k + 2e) + (Vw, b) 

+ AwV|,4| 2 (-2 + 3e) + Xu 2 \Vv\ 2 (a^ - 6) + c(e, A). 

Let us assume that k > is chosen so small that kv 2 < | in {u < 0}. This implies 
ip < 2v 2 . We may assume that A > 2u 2 in {u < 0} and get u 4 g < \Xu 2 ip\A\ 2 < 
Xu 2 v 2 \A\ 2 . We get 

4 ^ 6 ^ 6 < 
v 2 1 — kv 2 

Finally, fixing e > sufficiently small, we obtain 

(f t -A)w<(Vw,b)+c. 
Now, both claims follow from the maximum principle. □ 

Theorem 5.9 (C m+2 -estimates). Let u be as in Assumption 5.1. 
(i) There exists A > 0, depending on the C m+1 -estimates, such that 

tu 2 \V m A\ 2 + X\V m - 1 A\ 2 <c-X-t+ sup AlV^MI 2 

t = 
{u<0} 

at points where u < and < t < 1 . 
(ii) As in Theorem 5.8, initial smoothness is preserved. 

Remark 5.10. 

(i) This implies a priori estimates for arbitrary derivatives and any t > 0: It is 
known that estimates for u, v, \A\ and |V"M| 7 1 < m < M, imply (spatial) 
C M+2 -estimates for the function that represents the evolving hypersurface as 
a graph. Using the equation, we can bound time derivatives. 

(ii) For estimates at time t > I, we can use the previous theorems with t = 
replaced by t = to — 1/2. 
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(Hi) To control the m-th (spatial) derivative at time to > 0, we can apply the result 
iteratively and control the k-th derivatives, 1 < k < m, at time ^ . 

(iv) Theorem 5.9 implies smoothness for t > 0. We do not expect, however, that 
the decay rates obtained for | V"M| 2 are optimal near t = 0. 

Proof of Theorem 5. 9. Once again, we underline terms and factors that can be 
dropped to obtain uniform estimates up to t = 0. We define 

w :=tu 2 \V m A\ 2 + \\V m - 1 A\ 2 

for a constant A > to be fixed. We will assume that |V fe A| 2 is already controlled 
for any < k < m — 1. Suppose that < t < 1. The constant c is allowed to 
depend on quantities that we have already controlled. Thus the evolution equation 
for |V"M| 2 in Theorem 4.2 becomes for m > 1 

(f t - A) \\7 m A\ 2 < - 2 |V m+1 A| 2 + c \V m A\ 2 + c, 

(f t - A) \\7 m - l A\ 2 < - 2 |V"M| 2 + c. 



We get 



w = u 2 \V rn A\ 2 + 2t_uu \V m A\ 2 + tu 2 ^- \V m A\ 2 + A-^ jV™" 1 , 
Wl =2tuu t \V m A\ 2 +tu 2 (|V"M| 2 ) +A(|V I "- 1 A| 2 ) , 
= 2tu Uij \V m A\ 2 + tu 2 (\V rn A\ 2 ) + A ( |V m -M| 2 ) 

\ / ij \ / ij 

+ 2tu iUj \\7 m A\ 2 + 2tu ( Ui (|V m A| 2 ) + Uj (|V m A| 2 ) 



Wa = 



(| - A) w < u 2 \V m A\ 2 +tu 2 (-2 |V m+1 A| 2 + c \V m A\ 2 + cj 

+ A (-2 \V m A\ 2 + c) - 2t_\ Vu\ 2 \V m A\ 2 - ttu (vu, V | V rn A\ 
We observe that 

-ttu (Vu, V |V m ^| 2 \ <t_-\u\-c - \V m+1 A\ ■ \V m A\ < tu 2 I V m+1 A| 2 + c |V m A| 2 . 



Therefore we get 

( ; |~A) U ;<( C -2A)|V m A| 2 + C (A) 
and the result follows from the maximum principle. □ 

6. Holder estimates in time 

We will use the following Holder estimates to prove maximality of a limit of 
solutions. 

Lemma 6.1. Let u: M. n+1 x [0, oo) R be a graphical solution to mean curvature 
flow and M > 1 such that 

\Du(x,t)\ < M for all (x,t) where u(x,t) < 0. 

Fix any x € R n+1 and ti,t 2 > 0. If u(xo,ti) < —1 or u(x ,t 2 ) < —1, then 
\h-t 2 \ > 8{n+ \ )M 2 or 

\u(x , h )-u(x l t 2 )\ ^^Y ){M iy 
y/\h-t 2 \ 

The previous lemma implies that u is locally uniformly Holder continuous in 
time. Although Lemma 6 follows from the bounds for H provided by [11, Theorem 
3.1], we include below an independent and more elementary proof which employs 
spheres as barriers. 
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Proof. We may assume that t\ <t 2 . 

(i) Assume first that u(x ,ti) < —1. As \Du(x,t)\ < M for u(x,t) < 0, we 
deduce for any < r < 

u(x ,h) -Mr< u(y,t!) < u(x ,h) + Mr for all y e B? +1 (x ). 

Hence the sphere dB™ +2 (x , u(x , ii) + (M + l)r) lies above graph u(- 7 ti) and 
dB™ +2 {xQ, u(x , ti) — (M + l)r) lies below graph u(-, ti). When the spheres 
evolve by mean curvature flow, their radii are given by 

r(t) = y/r 2 -2(n + \){t-t{) 

2 

for ti < t < h+ 2 (n+i) ■ Both spheres are compact solutions to mean curvature 
flow. Hence they are barriers for graph u(-,t). In particular, we get 

u(x , h) - (M + l)r < u (x , h + 2 ^ +l ^ j < u{x , h) + (M + l)r. 

Set r := \[%n + l)(t 2 — ^i)- We may assume \ti — t 2 \ < 2 ( n +i)M 2 ■ Hence 
r < and the considerations above apply. We obtain 

u{x Q M) - (M+ 1)^2(71 + l)(i 2 - ti) <u{x Q ,t 2 ) 

<u(xo,h) + (M + l)y/2(n + l)(t 2 -h). 

Rearranging implies the Holder continuity claimed above. 

(ii) Assume now that u(xo,t 2 ) < —1 and u(xo,ti) > — 1. We argue by contradic- 
tion: Suppose that t 2 > t\ > t 2 — 8( - n+ 1 1 - )M 2 and 

(6.1) u(x , tl )-u(x ,t 2 ) ^ v ^ TI y (M + 1} 

V*2 — t\ 

Set r := y/2{n + l)(t 2 - t x ). We claim that 

(6.2) min{u(x , h), 0} - Mr > u(x , h) + r. 

If u(xo,ti) < 0, (6.2) follows by rearranging (6.1). Otherwise, we have that 
u(x ,t 2 ) + (M + l)r< - 1 + (M + 1)t/2(ti + l)(i 2 - *i) 



as M > 1. This proves claim (6.2). 

Now, using (6.2), we can proceed similarly as in (i): For some small e > 0, 
the sphere <9£>™ +2 (x , u(x , t 2 )+e) lies below graph u(-, t\) (for the positivity of 
e consider in (6.2) the terms —Mr near the center and +r near the boundary). 
Under mean curvature flow, the sphere shrinks to a point as t /• t 2 and 
stays below graphu(-,i). We obtain u(xo,t 2 ) +e < u(x n ,t 2 ), which is a 
contradiction. □ 

7. Compactness results 

Lemma 7.1. Let Q C B C R™ +2 and consider a function u: Vl — > R. Assume that 
for each a£R there exists r(a) > such that for each x £ fl with u(x) < a we have 
B r ( a ){ x ) flBcfi. Then £1 is relatively open in B and u{xk) — > oo if Xk — > x € dil, 
where d£l is the relative boundary of fl in B. 

Proof. It is clear that Q C B is relatively open. If u were not tending to infinity near 
the boundary, we find x n € f2 such that X r X €E dfl as n — > oo and u(x n ) < a for 
some a € R. Since £ r ( a )(x„) n B C fi, the triangle inequality implies x € B r ( a )(x n ) 
for n sufficiently large. This contradicts a; € dfl. □ 
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Remark 7.2. A continuous maximal graph is a closed set and - if sufficiently 
smooth - a complete manifold. 

Lemma 7.3 (Variation on the Theorem of Arzela-Ascoli). Let B C R n+2 and 
< a < 1. Let Ui : B — > R U {oo} for i G N. Suppose that there exist strictly 
decreasing functions r, — c: M. — > M + such that for each x G B and i > io(a) with 
Ui(x) < a < oo we have 

\uj{x) - Uj{y)\ 
\x - y\ a 



< c(a) for all y G B r ( a )(x) n B. 



Then there exists a function u: B — > M U {oo} such that a subsequence (u ik )ken 
converges to u locally uniformly in Q :— {x G B : u(x) < oo} and Ui k (x) — > oo for 
x G B \ O. Moreover, for each x G £1 wifft, u(x) < a we have -B r ( a +i)(x) flfic!! 
and 

|u(x) - 



F - 2/1 



< c(a + 1) for all y g B r ( a+1 ) (x) n -B. 



Proof. We adapt the proof of the Theorem of Arzela-Ascoli to our situation. Let 
L> := {x ( : / e N} be dense in B. 

If liminf Ui(xo) < oo, we choose a subsequence (ui k )keN, such that lim Ui k (x ) = 

i — ^oo k — >oo 

liminf m(xo). If liminf Uj(xo) = oo, we do not need to pass to a subsequence. 

Proceed similarly with X\, x 2l ■ ■ ■ instead of x - We denote the diagonal sequence 
of this sequence of subsequences by (Cti)ieN- Define u(xk) '■= lim Ui(xk) € MU{oo} 

i— ^oo 

for k e N. This limit exists by the construction of the subsequence (u»)ieN. By 
passing to the limit in the Holder estimate for Uj, we obtain the claimed Holder 
estimate with a + \ for u and a; = x k ,y — Xt, k,l £ N. Set u(x) := lim u(xk) 

k— >oo 

for x e B, Xk E D and — > x as k — > oo. The Holder estimate ensures that it 
is well-defined and fulfills the claimed Holder estimate with a + 1. Set il := {x € 
_B: u(x) < oo}. There, pointwise convergence and local Holder estimates imply 
locally uniform convergence in f2. □ 

Remark 7.4. 

(i) This result extends to families of locally equicontinuous functions. 

(ii) Notice that the functions Ui in the previous lemma are not necessarily finite 
on all of B. Hence the lemma can also be applied to functions Ui which are 
not defined in all of B: It suffices to set Ui := +oo outside its original domain 
of definition. 

(Hi) Observe that the domain f2 obtained in Lemma 7.3 may be empty. However, 
for the existence result (Theorem 8.2), the fact that f2 ^ is ensured by the 
choice of initial condition for the approximating solutions and Lemma 6.1. 



8. Existence 

In this section we will use approximate solutions to prove existence of a singu- 
larity resolving solution to mean curvature flow. 

We start by constructing a nice mollification of min{-,-}. Choose a smooth 
monotone approximation / of min{-, 0} such that /(x) = min{x, 0} for |x| > 1 and 
set min e {a, b} := ef (|(a — 6)) + b. 

We will set min e {w(x), L} := L at x if u is not defined at x. 

Lemma 8.1 (Existence of approximating solutions). Let A C R n+1 be an open set. 
Assume that u : A — > R is locally Lipschitz continuous and maximal. 
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Let L > 0, R > and 1 > e > 0. Then there exists a smooth solution R to 



u = VI + \Du\* ■ div ( J!^-= =] inB R (0) x [0,oo), 
\V 1 + l I)li l / 



u = L on &B fl (0) x [0,oo), 

k u(-,0) = min e {u a , e ,L} in B R (0), 

where u e is a standard mollification ofu . We always assume that R > Rq(L,s) 
is so large that L + 1 < uo,e 071 dB R (0). 

Proof. The initial value problem for R involves smooth data which fulfill the 
compatibility conditions of any order for this parabolic problem. Hence we obtain 
a smooth solution R for some positive time interval. According to [23], this 
solution exists for all positive times. □ 

Observe that the approximate solutions of Lemma 8.1 fulfill Assumption 5.1 with 

Cl={(x,t): u L e R < a} 
and there replaced by a for any a < L. 

Theorem 8.2 (Existence). Let A C R n+1 be an open set. Assume that u : A — > R 
is maximal and locally Lipschitz continuous. 

Then there exists fl C R n+1 x [0, oo) such that Ct n (R n+1 x {0}) = A x {0} and 
a (classical) singularity resolving solution u: fl — > R with initial value uq. 

Proof. Consider the approximate solutions R given by Lemma 8.1. The a priori 
estimates of Theorem 5.3 and Lemma 6.1 apply to this situation in {(x,t) £ B R (0)x 
[0, oo): u^ R (x,t) < L — l}. According to [11], we get uf^ i — > u L as i — > oo and 
m l is a solution to mean curvature flow with initial condition min{u, L}. 

Let us derive lower bounds for u L that will ensure maximality of the limit when 
L — > oo. As the initial value u fulfills the maximality condition for every r > we 
can find d — d(r) such that B r ((x, L — r — Vf) lies below graphminju, L} if |x| > d. 
Hence u L (x,t) > L — 2 for < t < ^jr — 2 (n+\) ^ \ x \ — d- Therefore for any 
T > there exists d > such that u L (x, t) > L - 2 for \x\ > d and < t < T. 

The estimates of Theorem 5.3, Theorem 5.8 and Theorem 5.9 survive the limiting 
process and continue to hold for u L : We get locally uniform estimates on arbitrary 
derivatives of u L in compact subsets of fin x (0, oo)) . The estimate of Lemma 

6.1 also survives the limiting process and we get uniform bounds for ||u L || co ,i ; o,i/2 
in compact subsets of SI. 

Now we apply Lemma 7.3 to u L , L E N, and get a solution (f2,u) and a subse- 
quence of u L , which we assume to be u L itself, such that u L — > u locally uniformly 
in Q. 

According to Lemma 7.1, Q is open in R n+1 x [0, oo). The C°' 1;0,1 / 2 -estimates 
imply that the domains of definition of uq and u\t=o coincide. In particular in 
Definition 2.1 we get A = ilo(Q) and u(-, 0) = u . 

The derivative estimates and local interpolation inequalities of the form 

||-EMIc°(B) ^ c(n,B) ■ \\w\\ C o(b) ■ \\w\\cz(B) 

for any w £ C 2 and any ball B (see e.g. [27, Lemma A. 5]) imply that u L — > u 
smoothly in n x (0, oo)). Hence u fulfills the differential equation for 

graphical mean curvature flow. 

The lower bound u L (x,t) > L — 2 above for \x\ > d and Lemma 7.1 imply 
maximality. 

Hence, we obtain the existence of a singularity resolving solution (f2, u) for each 
maximal Lipschitz continuous function uq : A — > M. □ 
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Remark 8.3. Notice that in the proof of Theorem 8.2 we started with the approx- 
imate solutions of Lemma 8.1 instead of u L in the proof of Theorem 8.2 as the 
former are smooth up to t = and allow to apply our a priori estimates. 

9. The level set flow and singularity resolving solutions 

In this section we explore the relation between level set solutions as defined at 
the beginning of Appendix A and singularity resolving solutions given by Theorem 
8.2. More precisely, we prove the following result 

Theorem 9.1. Let (Q,,u) be a solution to mean curvature flow as in Theorem 8.2. 
Let dT>t be the level set evolution of dflo as defined below. If &D t does not fatten, the 
measure theoretic boundaries of il t and T> t coincide for every t > 0: d^Ct t = d^T> t . 

For the definition of a level set solution and fattening, we refer to Appendix A. 

In order to prove Theorem 9.1 we need a few definitions which we summa- 
rize in Table 1. Unless stated otherwise, we will always assume that we consider 
signed distance functions which are truncated between —1 and 1, i.e. we consider 
max{— 1, min{d, 1}}, and negative inside the set or above the graph considered. 

(i) Let v: R n+1 x [0, oo) -> R be the solution to (A.l) such that v(-,0) is the 
distance function to dn a . Set V t := {x € R n+1 : v(x,t) < 0}. 

(ii) Let v. R n+2 x [0,oo) ->• R be the solution to (A.l) such that v(-,0) is the dis- 
tance function to 8Q x »■ Set C t := { (x, x n + 2 ) e R n+2 : v (x, x n+2 ,t) < 0}. 

(iii) Let w : W l+2 x [0, oo) — > R be the solution to (A.l) such that w(-, 0) is the dis- 
tance function to graphw(-, 0)|n o . Set E t :— {(x,x n+2 ) : w (x,x n+2 ,t) < 0}. 



solution to (A.l) 


initial set 


set 


w 


graph u Q 


Et 


V 






V 


on x R 


c t 



Table 1. Notation for weak solutions 



Theorem 9.1 will follow from 

Proposition 9.2. Let be a solution to mean curvature flow as in Theo- 

rem 8.2. If the level set evolution of <9Q does not fatten, we obtain H n+1 -almost 
everywhere that f2 t = V t for all t>0, i. e. V. n+1 (^t A2? t ) = for every t>0. 

We start by showing that v and v are closely related. 

Lemma 9.3. For v and v as above, we have v (x,x n+2 ,t) = v(x,t) for all points 
(x, x n+2 ,t) e R n+1 xRx [0, oo). This implies V t xR = C t and D+ x R = C t + . 

Proof. This follows directly from uniqueness of solutions to (A.l) as v (x, x n+2 , 0) = 
v(x,0). See Theorem A.l. □ 

Lemma 9.4. We have w > v. In particular, C C t + . 

Proof. This follows from w{-, 0) > v(-, 0) and Theorem A. 3. □ 
Lemma 9.5. We have graph u(-,t) C dE^ . 

Proof. Let w L : R n + 2 x [0, oo) -> R be the solution to (A.l) with w L (-, 0) equal to the 
distance function to graph u L , where u L is as in Theorem 8.2. According to [2] the 
solution w L does not fatten: For each e > there is a 5 > such that the inequal- 
ity w L (x, 0) > w L (x + ee n+2 , 0) + 5 holds if we truncate at appropriate heights. 
By Theorem A. 2 and Theorem A. 3 wc have that w L (x,t) > w L (x + ee n+2 ,t) + 
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S near the zero level set. Hence {(x,x n+2 ) e R n+1 x R : w L (x,x n+2 ,t) = 0} = 
graph u L (•, i). 

Observe that w L (-,0) w(-,Q). Hence Theorem A. 4 implies that w L (-,t) /- 
w(-,t) for all t > 0. 

Let x n+2 < u(x, t). Then x n+2 < u L (x, t) for some L and hence w L (x, x n+2 ,t) > 
0. Since w (x,x n+2 ,t) > w L (x,x n+2 ,t) > we have that 

(9.1) {(x,x n+2 ) : x n+2 <u{x,t)} C {(a;,a;" +2 ) : < w (x, x n+2 , t) } . 

On the other hand, for every (x, x n+2 , i) such that x n+2 = u(x, t) there is a sequence 
(x, it L (x, such that (a;, u L (x, t)) — > as L — > oo. Moreover, since the 

w; L converge monotonically, the convergence is locally uniform. We conclude that 

0= lim w L (x, u L (x, t), t) — w (x, u(x, t), t) . 

L— >oo 

This concludes the proof of graph u(-,t) C dE^ . 

By arguments similar to those used for proving (9.1), we can show that 

{(x,x n+2 ) : x n+2 >u(x,t)} C {(a;,^ 2 ) : w (x, x n+2 , t) < 0} . □ 
Corollary 9.6. Let x £ fit then w [x, x n+2 ,t) > for any x n+2 . 

Proof. The above argument in the case x n+2 < u(x,t) also extends to the case 
u(x, t) = +oo. □ 

Corollary 9.7. If C t or, equivalently, V t does not fatten, then fl t C T>^ . 

Proof. Combining Lemmata 9.3, 9.4 and 9.5 we get graph u(-,t) C T>f x R. This 
implies Q t C V+ . As V t is not fattening, we see that U n+1 (D+ \ V t ) = 0. Notice 
that V t C T>f C T>f. As fit is an open set, the claim follows. □ 

The following lemma shows that graph u(-, t) does not "detach" from the evolving 
cylinder at infinity. 

Lemma 9.8. We have T>t C fit- 
Proof. Denote by w R the solution to (A.l) with initial condition the distance func- 
tion to the set graph(w — R)- 

Notice that w R (-, 0) \ v(-, 0) as R — >• oo. Theorem A. 4 implies that 

(9.2) w R {-,t) \ v(-,t) as R -t oo. 

Suppose there are x, t such that x € T>t \ fit- Then by Corollary 9.6 it would 
hold for every R > and x n+2 that 

w R (x, x n+2 , t) > and v(x, x n+2 ,t) < 0. 

However, taking R — > oo this contradicts (9.2). □ 

Proof of Proposition 9.2. According to Corollary 9.7 and Lemma 9.8 we have 

VtCfl t C V? C V+. 

If there is no fattening U n+1 (X> t + \ V t ) = 0. The claim follows. □ 

Remark 9.9. 

(i) From Proposition 9.2 we have that 

sup {t > : u(-, t) ^ oo} = sup {t > : V t ^ 0} , 

i. e. the singularity resolving solution vanishes at the same time as the level 
set solution. Here u(x,t) — oo is to be understood as in Lemma 7.3. 
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(ii) Generically, level set solutions do not fatten, see [24]. Examples of initial 
conditions that do not fatten are mean convex hypersurfaces (see [30]) and 
star-shaped domains of definition (see [2] and references therein). 
(Hi) Under conditions similar to [2] it is possible to prove that w does not fatten 

and that (SI, u) is unique, 
(iv) Theorem 9.1 also holds if the 8SIq non-fattening assumption is replaced by 

non-fattening of the level set solution with initial condition graph uq. 
(v) If Dv along {v ^ 0}, we have — V t and hence Sl t —T> t . 

Appendix A. Definitions and known results for level set flow 

Different approaches have been considered in order to define a weak solution to 
mean curvature flow via a level set method (see for example [4, 12, 21, 28]). We 
define it as follows: Given an initial surface dEo, we define a level set solution to 
mean curvature flow as the set dE t = d{x : w(x,t) < 0}, where w satisfies in the 
viscosity sense the equation 

J ff - - WW) "V = in K"+ 2 x (0, oo), 
\«;(-,0) = to o (0 inR"+ 2 . 

and E = {x : w (x) < 0}. We also set E^ := {x: w(x,t) < 0}. 
We say that a solution to (A.l) does not fatten if 

H n+2 ({w(-,t) = 0}) = 

for all t > 0, where H n+2 denotes the (n + 2)-dimcnsional Hausdorff measure. 

Observe that our definition of solutions differs from the notion in [4, 12]: They 
define the level set solution to be {x: w(x,t) — 0}. If there is fattening, our 
definition picks the "inner boundary". Often, however, these definitions coincide, 
see e. g. [14, 21]. 

Let E C R n+2 be measurable. We define the open set E^, the measure theoretic 
interior of E, by 

E" := {x e K"+ 2 : 3r > 0: |B r (jc)| - \E n B r (x)\} . 

If E is open, we get E C E^ C E. We also define the measure theoretic boundary 
d^E of E by 

d"E := {x G K"+ 2 : Vr > 0: < \E D B r (x)\ < \B r (x)\) . 

In what follows we summarize some results in the literature that will be used 
in our proofs. We will work with the class BUC(Z) which are functions uniformly 
continuous and bounded in Z C E n+2 x [0,T). 

Theorem A.l (Existence [16, Theorem 4.3.5]). If w Q G BUC (E"+ 2 ) then there 
is a unique viscosity solution w G BUC (]R™ +2 x [0, oo)) to (A.l). 

Theorem A. 2 (Geometric Uniqueness [12, 16]). Let w\(x,t) and u>2(x,t) be vis- 
cosity solutions to (A.l) such that 

{x : Wi(x,Q) = 0} = {x : w 2 (x,0) = 0}, 

then 

{x : w\(x, t) — 0} = {x : w 2 (x, t) = 0} 

for any t > 0. 

Following Theorem 3.1.4 in [16] we have the following result for continuous sub- 
and super-solutions: 
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Theorem A. 3 (Comparison principle). Let w and v be continuous sub- and super- 
solutions of (A.l), respectively, in the viscosity sense in M. n+2 x [0,T). Assume 
that w and —v are bounded from above in R" +2 x [0,T). Assume that 

w(x,0) - v(x,0) < 

then 

t) - v(x, t) < for (x, t) e R n+2 x [0, T). 

Theorem A. 4 (Monotone Convergence [16, Lemma 4.2.11]). Consider functions 
fflo^m, wo € BUC(R n ) such that w , m /• w n . Then if w m and w are solutions 
to (A.l) with initial data w , m and w , respectively, we have for every time that 
w m /* w. 

Remark A. 5. 

(i) The (non-truncated) signed distance function to dE may be defined as dsix) = 
dist(x, E) — dist (x, R m \E). In particular, we assume that the signed distance 
function to dE is negative for every x e E. 

(ii) In general, the initial conditions considered in Section 9 will be given by trun- 
cated distance function to a set. 

(Hi) If the set dflo is compact and evolves smoothly under mean curvature flow, 
the level set formulation above agrees with the classical solution. 
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